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Abstract 

Brownian motion of a particle with an arbitrary shape is investigated theoretically. Analytical 
expressions for the time-dependent cross-correlations of the Brownian translational and rotational 
displacements are derived from the Smoluchowski equation. The role of the particle mobility center 
is determined and discussed. 


I. INTRODUCTION 


Brownian motion of particles with arbitrary 
shapes has been recently investigated in many 
different contexts, including proteins, DNA, 
nanofibers, actins or other biological nano and 
micro objects EN3- There is a rapidly growing 
number of experimental and numerical studies, 
which give rise to interesting fundamental ques¬ 
tions, challenging for theoreticians. 

For nanoparticles, the characteristic time tr of 
the rotational Brownian diffusion is typically much 
smaller than the time resolution t in the experi¬ 
ments, t » tr. Therefore, they can be treated as 
point-like spherical particles, and described by the 
standard Brownian theory [S|. 

However, for Brownian microparticles, tr is of 
the order of seconds, and therefore non-negligible 
in comparison to the typical time scales t of the 
measured Brownian motion, t < tr. For exam¬ 
ple, in Ref. 7J], the Brownian motion of a non- 
symmetric microparticle was investigated experi¬ 
mentally and numerically at the time scales com¬ 
parable with the characteristic time tr of the ro¬ 
tational Brownian diffusion, and in this case the 
standard approach [§] is not sufficient. Therefore, 
the idea was to measure and model numerically 
the time-dependent cross-correlations of the Brow¬ 
nian translational and orientational displacements 
of microparticles with different shapes. These re¬ 
sults were next also used to determine the mobility 
and friction matrices. 

The goal of our work is to derive theoreti¬ 
cally the explicit analytical expressions for cross¬ 
correlations of the Brownian translational and ori¬ 
entational displacements at all the time scales, and 
for arbitrarily shaped particles. 
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Parts of such a theoretical analysis have been al¬ 
ready done. The motivation came from Ref. Q, 
focused on the problem of determining whether 
measurements of the self-diffusion coefficient and 
intrinsic viscosity of fibrynogen can be used to de¬ 
termine the protein configuration. Evaluation of 
both these quantities for particles with arbitrary 
shapes hinges on solving fundamental problems: 
what is the Brownian contribution to the intrinsic 
viscosity, and is the self-diffusion coefficient sensi¬ 
tive to the choice of a point on particle for which 
the mean square displacement is determined. 

The expressions for the intrinsic viscosity of a 
Brownian non-symmetric particle were theoreti¬ 
cally derived in Refs. id- The relation between 
the Brownian mean square displacements of differ¬ 
ent points of a particle was determined in Ref. [llj. 

In this work, we consider a single Brownian par¬ 
ticle of an arbitrary shape, in general non-isotropic 
and non-axisymmetric. Starting from the Smolu¬ 
chowski equation [12]- [l4|, we develop a new for¬ 
malism, which allows to determine the particle ro¬ 
tational and translational motion in a much sim¬ 
pler way as e. g. in Ref. [lfj, which is based on 
the Euler angles and Wigner functions. 

Moreover, the essential result of this work is 
that using our new formalism, we derive sim¬ 
ple explicit analytical expressions for the cross¬ 
correlations of the Brownian translational and ro¬ 
tational displacements. No such formulas have 
been known yet - instead, numerical Brownian sim¬ 
ulations have been extensively used, as e. g. in 
Ref. 0- 


II. SYSTEM AND ITS THEORETICAL 
DESCRIPTION 

We consider an isolated Brownian particle of ar¬ 
bitrary shape immersed in an unbounded fluid of 
viscosity rj and temperature T. Its state will be 
described by X = (R, O), where R is the position 
of the chosen particle point and fl describes the 
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particle orientation, in particular in terms of the 
Euler angles which specify the orientation of the 
particle body-fixed axes with respect to the space- 
fixed axes. Probability distribution of finding the 
particle at time t in state X will be denoted by 
P(X,t). It is normalized as 

JdXP(X,t) = 1, (1) 

with 

dx = dn dh. ( 2 ) 

and the element dfi defined in terms of the Euler 
angles in the same way as in Ref. on page 161. 

The evolution of the probability distribution 
P(X, t) is governed by the Smoluchowski equation 

^P(X,t)=C(X)P(X,t) (3) 

with the Smoluchowski operator, which in the ab¬ 
sence of external fields has the form 

C{X)= Vx -D(h)-v x . (4) 

In the above equation, 

< 5 > 

with 


d 

f d d d \ 

(6) 

m 

V<TRi ’ dlh ' dB-J ’ 

d 

( d d a \ 

(7) 

da 

\c>a:i ’ da .2 ’ da .3 ) 


where ctk is the angle of rotation around the axis fc; 
here and in the following, we denote the Cartesian 
indexes corresponding to x, y, z by k = 1, 2, 3. Ac¬ 
cording to the fluctuation-dissipation theorem 0], 
the 6x6 diffusion matrix D is proportional to the 
6x6 mobility matrix /x, 

£>(«) =k B Tn{n), ( 8 ) 

It is important to stress that the matrix /x does 
not depend on position R, and its dependence on 
follows entirely from the particle rotation; it 
means that /x rotates with the particle (see Ap¬ 
pendix IB]) . As we see from Eq. (|8]l. the single¬ 
particle Brownian motion is determined by its 6x6 
mobility matrix /x. If applied to the external force 
F and torque T exerted on the particle, in the ab¬ 
sence of external ambient flows, it results in the 
particle translational and rotational velocities, U 
and uj, 



The mobility matrix fi consists of four Cartesian 
3x3 sub-matrices, 

with the indices t, r denoting the translational and 
rotational components, respectively. The 6x6 mo¬ 
bility matrix fi is symmetric [IH, therefore /x-j = 
= tfji- Therefore, the 3x3 
matrices fi u and /x rr are symmetric, but /x tr in 
general is not symmetric. Analogical notation is 
adopted and the symmetry properties are held for 
the corresponding 3x3 sub-diffusion matrices, de¬ 
fined with the use of Eq. (0). 

In general, the mobility matrix /x depends on the 
choice of a reference center which is observed. The 
hydrodynamic mobility center is such a point for 
which the rotational-translational mobility matrix 
H rt is symmetric. The position of this point is 
explicitly specified e.g. in Refs. 0 and list . 

Transformation relations between the mobility 
Hi for the reference center iq and the mobility /x 2 
for the reference center r 2 , called translational the¬ 
orems for mobility matrices, are the following, 

H r 2 r = Hi = M rr , (11) 

M? = + M rr x(r 2 -n), (12) 

H 2 = M? - (r 2 - rr) x x (r 2 - rq) 

- (r 2 - ri) x H r i + Vi x ( r 2 - ri), (13) 

with the notation that for a second rank tensor A 
and a vector V, one has 

[AxV ]ij = AiktjkiVi, [V x A]ij=eukViAkj. (14) 

Here and later on, we use the Einstein’s summation 
convention, unless it is explicitly written that the 
repeating indices are fixed. 

In particular, the translational theorems listed 
above can be applied to obtain the mobility ma¬ 
trices for an arbitrary reference center, if He for 
the hydrodynamic mobility center is known. The 
advantage of using this special center is not only 
the simplicity of the corresponding expressions. In 
Ref. [IJJ it has been shown that the hydrodynamic 
mobility center is especially meaningful for the 
translational diffusion process. In this work, we 
will show that such a reference center is also es¬ 
sential for the analysis of the cross-correlations of 
the Brownian translational or rotational displace¬ 
ments. 

Taking into account the structure of the mobility 
and diffusion matrices, we split the Smoluchowski 
operator, defined in Eq. into four parts, 

C = C tl + L tr + C rt + C rr 


(15) 
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with 


C tt = 

9 D u d 
dR k kJ dR 3 ' 

(16) 

£ rr = 

d d 

da k kj d aj 

(17) 

c rt = 

d Drt d 
da k kj dRj ’ 

(18) 

C tr = 

5 D tr 9 

8 R k k3 da 3 ' 

(19) 


In the above expressions for C u and C rr , the order 
of the diffusion coefficients and both derivatives 
can be interchanged. This property can be easily 
derived taking into account that D u does not de¬ 
pend on R, and D rr is symmetric. Below, we will 
show that for any second rank symmetric tensor S, 

A. S = °. ,20) 

For the choice of the hydrodynamic mobility cen¬ 
ter, D rt and D tr are also symmetric and therefore 
in Eqs. m-m the order of all the operations is 
arbitrary. 

We will now derive the relation (0(71) . Let B be 
a tensor of rank n. Suppose that B rotates with 
the particle, and the rotation matrix TZ transforms 
an initial particle orientation $2o to a final orienta¬ 
tion Cl. Then, the components -Bi 1 t 2 ...i„ of B, with 
i k = 1,2,3 for k = 1, ...n, change according to the 
formula 

-®*i»2...i» (^) = 'R'iiji'R-iaj? (^o)i (21) 


where 1Z; are the components of the rotation ma¬ 
trix n. 

One can show Q that 


da k 1112 


— e kj 1 i 1 Bj 1 i 2 . 


' e kj 2 i2^ilj2- 


( 22 ) 


It follows from the above that for a second rank 
symmetric tensor S the r.h.s. of Eq. (1221) is equal 
to zero, what means that Eq. m is satisfied. 


III. BROWNIAN TRANSLATIONAL AND 
ROTATIONAL DISPLACEMENTS 


In this section, we will describe Brownian mo¬ 
tion of a particle with an arbitrary shape. At time 
(=0 we choose an arbitrary particle point R(0) and 
orientation $1(0), and we trace the positions R(f) 
and orientations Cl(t) at times t. 

We first analyze the time-dependent Brownian 
translational displacements, 


(23) 


We want to evaluate the dynamical cross¬ 
correlations, 


(AR(f)AR(f)) 0 = 

I dX (R-Ro)(R-Ro)P(R,n 


Rr 


,n 0 ;t) • (24) 


The average (...)o is taken with respect to the par¬ 
ticle positions R and orientations fl at time t, with 


the conditional probability P^R, fl H a ,Cl 0 ;tj, 
which takes into account that at t = 0 the par¬ 
ticle is located at R(0) = Ro, and oriented along 
$1(0) = Cl Q . Actually, the average depends on R 0 
and Cl 0 and formally, it should be explicitly de¬ 
noted as (...) R q ) but we use a simpler notation 
(...)o- The average is taken with respect to the 
conditional probability, which satisfies the Smolu- 
chowski equation ([3j), i.e. it has the form 


p(r,« 


Ro,$lo;t^ = exp(£f)d(R — R 0 )<5($1—fi 0 )- 


(25) 

To determine change of the particle orientation 
Cl, we follow Ref. 0 and at time t = 0 we in¬ 
troduce three mutually perpendicular unit vectors 
u*74 (0) = Ug p ' ) which characterize the particle ori¬ 
entation, with p = 1,2,3, and we trace their evo¬ 
lution in time. 

The orientation u*P (f), p = 1, 2, 3, at time t can 
be interpreted as the result of a rotation matrix 
7Z(t), acting on the initial orientation ub4(0), 


U W(f)=7l(^uW(0). (26) 

With the use of the initial basis 11*74(0), the com¬ 
ponents of the rotation matrix are expressed as 
IZpq = u( p )(0) • u ( q \t), with p, q = 1, 2, 3. 

The matrix 7Z is orthogonal, i.e. 1Z T = "R -1 , 
where the superscript T stands for the transpo¬ 
sition. We now decompose the rotation matrix 
into symmetric, 7ZS s \ and antisymmetric (skew- 
symmetric), 7Z^ a \ parts, see Eqs. (IA4I) - (IA5I) . 


Ti = n is) + n {a) . ( 27 ) 


The antisymmetric part can be used to construct 
a vector, defined as 

A u k (t) = (28) 

This vector is parallel to the Euler rotation axis 
0, with its length equal to the absolute value of 
the sinus of the rotation angle around this axis. 
From Eq. (l26l) it follows that 


Au(f) = ^uW(0)xuW(f). (29) 

P= 1 

Following Ref. Q, in this paper we use Au(f) 
to describe the rotational Brownian motion. It is 


AR = R(t) — R(0). 
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good to keep in mind that the antisymmetric part 
n (°) 

is not sufficient to describe all the properties 
of the rotation E3- 

The conditional probability of the occurrence of 
a given rotation has the form, 


p(ft ft 0 ;tj = j ft Ro, flo; tj . (30) 


It satisfies the Smoluchowski equation with the ini¬ 
tial condition P(J~l, t = 0) = 8(ft — flo). Therefore, 


p(n n 0 ;t) =exp (C rr t)S(fl 


fin 


(31) 


The conditional probability depends only 
on the rotation matrix 7Z. which leads from fl 0 to 
fl (it does not depend on specific initial orienta¬ 
tion fl 0 ). The Smoluchowski operator C rr is self- 
adjoint [8|], and the forward and backward Smolu- 
chowski equations are the same. Therefore, the 
probability associated with 7Z is equal to the prob¬ 
ability associated with TZ~ l . Saying it differently, 
this property follows from the detailed balance con¬ 
dition. [s] 

Taking this property into account, we conclude 
that the symmetric and antisymmetric parts of the 
rotation matrix are not correlated, 


('R.^ s) (t)7l {a) (t)) 0 = 0. (32) 

To derive Eq. (1571) . we decompose 7 Z~ x = 1Z T into 
symmetric and antisymmetric parts, 

-R.- 1 =1Z {S) -lZ {a \ (33) 

and use the relations 

(n^(t)n^(t)) 0 = ((n- 1 )^(t) (iz- 1 )^(t)) 0 
= -(n {s \t)n {a \t))o (34) 

In Ref. Q, the Brownian motion was described 
based on measurements of the time-dependent 6 x 
6 cross-correlation matrix, 


C(t ) 


(AR(f)AR(f)) 0 (AR(f)Au(f)) 0 
(Au(f)AR(f)) 0 (Au(t) Au(f)) 0 


■ (35) 


The diffusion matrix D was determined as the 
time-derivative of the correlation matrix C (t) at 
time t = 0, 


1 

2 



t—0 


= D, 


(36) 


with the components, in analogy to Eq. (1TT)1) . de¬ 
noted as, 


D u D tr 

j~) r t j^rr 


(37) 


The relation Eq. (157)1) follows from the Smolu¬ 
chowski equation ©. 


In this work, we will describe Brownian mo¬ 
tion of a particle with an arbitrary shape by de¬ 
riving explicit analytical expressions for the cross¬ 
correlation matrix C (t) at all times t. In general, 
C (t) depends on the choice of the reference cen¬ 
ter. The transformation relations for Eq. (135|) to 
another reference center depend on both the an¬ 
tisymmetric and symmetric parts of the rotation 
matrix, 7 Z^ and TZ ^, respectively. These rela¬ 
tions will be discussed in Appendix [Cl 


IV. EVOLUTION DUE TO ROTATIONAL 
DIFFUSION 


A. Tensors which rotate with the particle 


While evaluating time-dependent cross¬ 
correlations, it is important to know how to 
determine changes due to the rotational diffusion. 
In this case, it is sufficient to average with respect 
to P(ft\fl 0 ;t), given by Eq. (15T1) . 

Therefore, we consider a tensor B which rotates 
with the particle and assume that ft(t = 0) = flo. 
For the conditional average of this tensor we find 


(B(£)) 0 = J dfl B(fl) P(ft fl 0 ;tj 

= J dft B(ft)exp(C rr t)8(ft-ft 0 ) 

= f dft {exp(r rr f)B(0)}ci(n-Oo) 


= exp(£ rT, f)B(fi) 


£2—fin 


(38) 


The explicit form of the above average depends 
on the tensor rank. It can be relatively easy eval¬ 
uated in the body-fixed frame in which the rota¬ 
tional diffusion tensor D rr is diagonal, 


DP = DiSij , for given i,j = 1, 2, 3. (39) 

From now on, this frame will be used until the end 
of the paper. 

The time-dependence of a tensor B follows from 
solving the eigenproblem for the Smoluchowski op¬ 
erator C rr , defined in Eq. (1171) . The eigenval¬ 
ues depend on the rank of this tensor. For ex¬ 
ample, when the tensor rank is one (i.e. when B 

is a vector), we find three eigenvalues — f- 1 ' 1 (as in 
Ref. [H|), with 

/P = 3D-D U t = 1,2,3, (40) 

where 

D=l(D 1+ D 2 + D 3 ). (41) 

Keep in mind that in this work, and D always 
refer to the rotational-rotational diffusion (we skip 
the superscript rr for simplicity). 
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The corresponding eigenvectors are parallel to 
the principal axes of the operator D rr . Denoting 
by T)(0) the components of this vector in the frame 
of reference in which D rr is diagonal at t = 0, and 
taking into account Eq. (l38l> . we obtain, 

W(t))o = exp(-if»t)VS(0). (42) 

For a second rank tensor H(t), the problem was 
solved in Ref. and the results are listed in Ap¬ 
pendix lAl The antisymmetric part of H{t) may be 
associated with a vector, and therefore its evolu¬ 
tion follows from Eq. (PUZll . with the same eigen¬ 
value given by Eq. flUl) . The evolution 

of the traceless symmetric part of H(t) is associ¬ 
ated with five other eigenvalues, — f!?\ and 

—as the characteristic exponents in the time 
decay, where 


jf } = 3 {Di + D), i = 1,2,3, (43) 
/(+) = 6D + 2A, (44) 

/(-) = 6.D — 2A, (45) 


with Di and D defined in Eqs. m and (ED and 

A = yjD\ + D\ + — D ± D 2 - D!D 3 - D 2 D 3 . 

(46) 

The trace of H(t) is a scalar and as such, it does 
not depend on time. 

The expressions for the time-dependence of the 
average ( H(t))o , analogical to Eq. (l42l) . are listed 
in Appendix [A] in Eqs. (IA2I) - (IA3I) . The important 
outcome of our analysis is that these expressions 
contain different characteristic exponents, listed in 
Eqs. (TO and fi3l) - (H5l) . 


B. Time scales of the rotational diffusion 


From our analysis, outlined in the previous sub¬ 
section, it follows that in general, there is no single 
characteristic time scale tr of the rotational self¬ 
diffusion. The exponents, listed in Eqs. pop and 
m (USD, determine several characteristic time 
scales of the translational and rotational correla¬ 
tions, 


V. MEAN PARTICLE DISPLACEMENT 


In general, due to the rotational diffusion, the 
mean particle displacement is not equal to zero i.e. 
if R is the position of the chosen particle point 
(AR(f)) 0 = ((R(t) - R(0))) 0 / 0. To find this 
mean displacement let’s first consider special case 
when the particle point is the mobility center Rc- 
Let’s calculate the time derivative of (ARc(f))o- 
Taking into account the Smoluchowski equation 
© we have 


^< AR c(% = 
f dX (Rc~Rco)£-P (l 


R c , J~I 


R-co, ^o; t ). (48) 


Now we will use the property, that for the mo¬ 
bility reference center, all the derivative can be 
shifted left. Then, integrations by parts give that 
the above derivative is equal to zero. Thus, 


(AR c (f)) 0 = 0. (49) 


Now consider the mean displacement for arbi¬ 
trary chosen particle point R. The vector R(t) — 
Rc(t) rotates with the particle. Thus from Eq. 
m we obtain 


((R(t) - Rc(t))) 0ii = (R(0)-Rc(0))i exp(-/( 1) t). 

(50) 


and with (TTO we get 


(AR(l)) 0 ,i = (Rc(0) ^ R(0))i {l - exp(-/«f)} 

(51) 

From Eq. (15T1) it follows that for long times, the 
mean position of an arbitrary point of the particle 
tends to the same limit: the initial position of the 
mobility center, 


(R(t)) 0 -A Rc(0), for t —> oo. (52) 

This result can be used to experimentally deter¬ 
mine the particle mobility center, what is espe¬ 
cially useful if the the location of i?c(0) is not 
known in an analytical form. 

In sections IVI1IVIIII we will follow the motion 
of the mobility center, R(t) = Rc(f), chosen as 
the reference center in the mobility matrix, see 

Eqs. CEB-CES. 


= 1 // (+) , T ( ) = l// ( 

= 1/jf \ rf : > = l//i 2) , (47) 


VI. TRANSLATIONAL-TRANSLATIONAL 
CORRELATIONS 


with i = 1,2,3. In general, these scales differ 
from each other, and a careful analysis of the time- 
dependence at long times in needed. 

It is always true that t - 1 ' 1 >t- 2 ' > and >t^ + \ 
but can be smaller or larger than t^~\ de¬ 
pending on the particle geometry. 


A. General expressions 

Let us now consider correlations of a Brownian 
particle displacements, (ARc(t)ARc(i)) 0 . We re¬ 
mind that we stay in the reference frame in which 
D rr is diagonal and follow the motion of the mo¬ 
bility center Rc. 
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It is convenient to calculate first the time deriva¬ 
tive of the expression (EH) : 


dt 


(ARc(f)ARc(f)) 0 = 


fdX (R c -R co )(R c -R co ) C P (R c , ^ Rco, ^o! t) 


(53) 


Taking into account the explicit expressions 
m-m for the Smoluchowski operator and in¬ 
terchanging the order of derivatives we obtain, 

j t (AR c (t)AR c (f)) 0 = 2 (D tt (t)) 0 . (54) 

We perform the above average using Eq. 

Since D u is a symmetric tensor we get from Eq. 
(IA2I) the following result for diagonal components, 

lj t (^c(t)A-Rc(t)) 0 ,« = 

D* + e- /(_>t I?“ (_) + e~ f(+h Dll (+) (55) 


square displacement is a scalar, and therefore it 
does not depend on the orientation, and in partic¬ 
ular, on the initial orientation fio- 

The mean square displacement of the mobility 
center, given in Eq. reproduces the result 

from Ref. HU- 

For long times, Eqs. m and (1591) take the form, 

\ <ARc(i)ARc(t)> 0 ,« = DH + 

+ + (61) 

\ (ARc(t)ARc(t)) 0 , y = 4y D% + 0{e-^). 

Jk 

(62) 

It is important to emphasize that there are con¬ 
stant non-vanishing terms in the above expres¬ 
sions. 


B. Special cases 


where z = 1,2,3 and 

D* = -tr D tt . 
3 


For an axisymmetric particle, a frame is chosen 
where the rotational-rotational diffusion matrix is 
(56) diagonal, with the coefficients 


The meaning of the other symbols is explained in 
Appendix E] After integration of Eq. (l55ll with re¬ 
spect to t one obtains for i = 1,2,3, 


\ (ARc(t)ARc(t)) 0 ,« 


l-e-/ (+,t 


/(+) 


T-ytt ( + ) 


dH + 


l-e-/ ()t 


D 


tt (-) 


(57) 


D X = D 2 ^ D 3 . (63) 

The translational-translational diffusion matrix is 
also diagonal, 

Dfj = D\5ij, for i,j = 1,2,3, (64) 

D\=D\± D\. (65) 

The diagonal correlations have the form, 


The off-diagonal components i ^ j are given by 
~ (AR c (f)AR c (f)) 0 .. = e-f^DH (1)8) 

where k is the remaining third index, k ^ i,j. Af¬ 
ter integration of Eq. (1551) with respect to t one 
obtains 

1 i _ e -/( 2) t 

- (ARc(t)ARc {t))o,ij = (2) D%. (59) 

Jk 

In the limit of t —> 0, the time derivative of the 
expressions (EZD-HMD approaches the correspond¬ 
ing elements of the diffusion matrix, in agreement 

with Eq. (1551) . 

From the relation (1571) , taking into account that 
D n and D u ) are traceless, one derives the 
expression for the mean square displacement, 

\ (ARc(f) • ARc(t)) = D t t. (60) 
6 

In Eq. (1501) . the subscript “0” associated with the 
averaging has been omitted, because the mean 


2 (ARc(t)ARc(t)) 0iij = 

1 _ P —6Dit 

Dt t+ ^ (Dl-D*), (66) 

for *=1,2, and 

2 (ARc(t)ARcr(t)) 0i 33 = 

1 _ p—6Dit 

D f t -(67) 

The off-diagonal components vanish. 

For a spherical particle, the rotational and trans¬ 
lational diffusion tensors are isotropic, 

D x = D 2 = D 3 = D, (68) 

D\ = D\ = D\ = D l , (69) 

and 

^ (ARc(f)ARc(f)) 0 = D f t I, (70) 

where I is the identity tensor. 
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VII. ROTATIONAL-ROTATIONAL 
CORRELATIONS 

A. General expressions 

From Eq. (EH1) it follows that the rotational- 
rotational correlations read, 

(Au(t)Au(t)) 0 = 

-lEE u o P) x (u {p \t)u^(t)) o x u^. (71) 

P= 1 9=1 


For a spherical particle, Di = D 2 = D 3 = D, and 
(Au(t)Au(t)) 0 = 


I _ lp-6 Dt , 1-2 Dt 
6 12 6 + 4 6 


I. 

(77) 


VIII. ROTATIONAL-TRANSLATIONAL 
CORRELATIONS 


A. General expressions 


The average (...)o in the above equation is eval¬ 
uated based on the expressions provided in Ap¬ 
pendix [A] 

The off-diagonal elements vanish, and the diag¬ 
onal elements, * = 1,2,3, are given by 


(Au(f)Au(t)) 0|ji = ^ - 

3(A D) + A _^(+) t 

12A 6 


3(.D — Di)-\- A 

- -f J 

12A 

1 A2) 1 f (l) 

__ e -/, 4 + -e-^ 4 . (72) 


In the limit of t —>• 0, the time derivative of the 
above expressions approaches the corresponding el¬ 
ements of the diffusion matrix, in agreement with 
Eq. m- 

For long times, 

lim (Au(t)Au(t)) 0 = ^1. (73) 

t—>-oo o 

This result can be also obtained as the average of 
Eq. (1721) with the equilibrium distribution. 


Using Eq. (l29l) we have, 


(Au(f)AR c (t)) 0 = i Y,u (p) x ^u w (t)AR c (i)^. 

P= 1 

(78) 


where, as before, Uy^ with p = 1,2,3 are the unit 
vectors corresponding to the orientation L2o, and 
we choose the body-fixed frame in which the rota¬ 
tional diffusion tensor D rr is diagonal. 

As shown in Appendix (Bj 

<J U ( p) (t)AR c (t)) o = -2 J dTe- f p 1){t ~ T) (A>)(r)^, 

(79) 


where the second rank tensor A^ is given by 


a (p) = u (p) x D rt . (80) 


B. Special cases 


For an axisymmetric particle, Di=D 2 ^D 3 , and 


(Au(t)Au(i)) 0 n = (Au(t)Au(t)) 0 22 = 


/O.ll ~~ \^- 1Ll V < y^- 1L H l y/o,22 
1 1 . 1 1 
6 


1 _ lg-6Dit _ £ e -(5Di+D 3 )t _|_ _ e ~(D 1 +D 3 )t 


6 4 

(Au(f)Au(f)) 

1 


4 


—2D\t 


— 1 i _ 1 -(2D 1 +4D 3 )t 

°’ 33 “ 6 + 12 6 2 6 

(75) 


To describe the rotational self-diffusion of an ax¬ 
isymmetric particle, it is common [Til to trace the 
change of the particle orientation vector along the 
symmetry axis, in this paper denoted as u (, ^(t). 
This vector rotates with the particle, and there¬ 
fore its evolution follows from Eq. (1421) . Using this 
relation, we reproduce the standard formula 1^], 




= 2(1 - e~ 2Dlt ). 


(76) 


We now decompose the tensor A( p ) and deter¬ 
mine (A( p )(f)) 0 for all its parts as described in Ap¬ 
pendix m The resulting expressions are plugged 
into Eq. G9, and the integration with respect to 
r is performed. 

We obtain the following diagonal Cartesian 
components of the rotational-translational corre¬ 
lations, 


(Au(t)AR c (f)) 0ii = - 


D lh D ll 


A 


8 Dk 

D%+D% + _ /f D t 
-— te Ik 


D$+D 


kk 


t e 


fPt 


(81) 


where i = 1, 2, 3 and j, k are the remaining second 
and third indices such that j + k = 6 — *. 

The off-diagonal Cartesian components are, 














/O ,ij 


-f 

e —e 3 

~7F z fF 


(!) t 


+ 


D rt 

ij 

Dj-Dk + A e~f ( h 

2 

l 

> 

T 


— e 




A i) 
Ik 




Dk — Dj + A e - /t+>t -e-^ 


/{+) _ / ( 1) 


* ± 3, (82) 


where the remaining third index k = 6 — i — j. 
The expression (1821) is not symmetric in ij. but its 
derivative at t = 0 is symmetric. 

In the limit of t —> 0, the time derivative of the 
above expressions approaches the corresponding el¬ 
ements of the diffusion matrix, in agreement with 
Eq. (1M1) . For long times, 

lim (Au(t)ARc(t)) = 0. (83) 


B. Special cases 

For an axisymmetric particle, D tr is skew- 
symmetric jTI[ . Therefore, it vanishes if evaluated 
with respect to the mobility center. Taking this 
into account, we obtain from Eqs. (I5T1) - (1521) the 
simple result, 

(Au(f)AR c (t)) 0 = 0 . (84) 

IX. CONCLUSIONS 

In this work, we performed theoretical analysis 
of the Brownian motion of a particle with an arbi¬ 
trary shape. We derived analytical expressions for 
the time-dependent cross-correlations of the dis¬ 
placements of the particle position and orientation. 
These results were written in the frame of reference 
in which the rotational-rotational diffusion tensor 
is diagonal at t = 0, and with the choice of the mo¬ 
bility center as the reference center. These results 
can be compared with experimental data, using the 
following procedure. 

Based on measurements of the time-dependent 
orientation u^^(t), it is possible to determine the 
correlation tensor (Au(t)Au(f)) 0 and find its prin¬ 
cipal axes. From our analysis it follows that they 
coincide with the principal axes of the rotational- 
rotational diffusion tensor. Therefore, to be com¬ 
pared with our theoretical expressions, the exper¬ 
imental data should be recalculated to this new 
frame of reference in which (Au(t)Au(f))o is diag¬ 
onal. 

The initial position of the mobility center Rc (0) 
can be determined experimentally by tracing in 
time the average position (R(t))o of any particle 
point, and using our finding that (R(f))o —5- Rc(0) 
when t —>■ oo. 


The transformation formulae for C(f) from 
the mobility center to another arbitrary point 
can be derived analytically, as described in Ap- 
pendix[Cj However, the transformation formula for 
(AR(f)AR(f))o involves the symmetric part of the 
rotation matrix TiS s \ which cannot be expressed 
in terms of C (t). 

Our simple explicit analytical expressions for the 
correlations C (t), valid for the mobility reference 
center, can be compared with the experimental 
data, measured for an arbitrary reference center, 
using the following procedure. Once the princi¬ 
pal axes of the rotational-rotational diffusion ten¬ 
sor are determined, and used as the new coordinate 
system, with the recalculated time-dependent ori¬ 
entation ( t ), and the initial position of the mo¬ 

bility center Rc(0) is determined experimentally, 
this information can be used as follows. First, the 
difference between the initial position R(0) of an 
arbitrary reference center, traced in experiments, 
and R c (0), can be expressed as a linear combina¬ 
tion of the orientation vectors u( p )(0), 

3 

R(0)-R c (0) = ^a p uW(0), (85) 

p=i 

and the coefficients a p can be determined. Tak¬ 
ing into account that R(f) — Rc(i) rotates with 
the particle, we can use the same coefficients a p to 
evaluate Rp(t) as 

3 

R c(t) = R (t) - ^2 a p u (p \t), (86) 

p= i 

Eq. (1861) allows to express (the unknown) 
stochastic trajectory Rc(t) in terms of (the 
known) stochastic trajectory R(f) and orientation 
u^(t). This allows to extract from the mea¬ 
sured data the correlations of the Brownian dis¬ 
placements of the mobility center, in the frame in 
which the rotational-rotational diffusion matrix is 
diagonal, and compare with our theoretical expres¬ 
sions for the time-dependent cross-correlation ma¬ 
trix C(f) for the mobility center Rc(t). 

In contrast to numerical simulations, the ana¬ 
lytical expressions provided in this work allow to 
determine the cross-correlations exactly. The ac¬ 
curacy of theoretical expressions is especially im¬ 
portant for times comparable to the characteristic 
time scales of the rotational diffusion, when the 
cross-correlations change significantly with time. 











9 


Acknowledgments 

M.L.E.-J. and E.W. were supported in part by 
the Polish National Science Centre under Grant 
No. 2012/05/B/ST8/03010. M.L.E.-J. benefited 
from the scientific activities of the COST Ac¬ 
tion MP1305. 


Appendix A: Change of a second rank tensor 
due to the rotational diffusion 

In Eqs. (I42l) - (l40ll . we derived the time- 
dependence of a vector, which rotates with the 
particle. Now we will do the same for a second 
rank tensor. 

Let H(t) be any second rank tensor which ro¬ 
tates with the particle, with H = H( 0). From Eq. 
(1551) it follows that time evolution of this tensor 
due to rotational diffusion only is given as 

(H(t)) o = exp (C rr t)H, (Al) 

The explicit expression for H(t) follows from Eqs. 
m (see Ref. @). As in the whole paper, we adopt 
the frame of reference in which the rotational- 
rotational diffusion tensor is diagonal, as in Eq. 
®. For the diagonal components we obtain, 

(• H(t))o,u = ^tr H + exp(—/ (+) t)JT[ +) 

+ (A2) 

where the first term is a scalar, and as such, it 
does not change in time, and the characteristic ex¬ 
ponents /1 + ) and /are given in Eqs. HU' 03' 


j 


For the off-diagonal components i ^ j, 

i H (t))o,ij = exp(-fl 1] t)Hlf + exp(-/f : } t)H[f 

(A3) 


where k is the remaining third index, k ^ and 
the characteristic exponents f- * ^ and are given 
in Eqs. (I40l)-(l43l). 

In Eq. (IA3I) , we use the decomposition of a sec¬ 
ond rank tensor H into symmetric and antisym¬ 
metric parts, defined as 



H {a) 


ij 


2 (Hij + Hji) j 
2 (-®»J — Hji) ■ 


(A4) 

(AS) 


In Eq. (IA2I) , the symmetric diagonal part is fur¬ 
ther split into the isotropic part, |tr H 5ij, and the 
diagonal non-isotropic parts, 




a 


'1 3 Dt-D\ 

2 T 4 A ) 



± Dj 4A Dk (Hjj-Hkk). (A6) 
Above, j and k ^ j are the remaining indices 
different than i. 


Appendix B: Rotational-translational 
correlations 

Using Eq. (BSD, one can write 


(uW(f)AR c (f)) o 



dRc (Rc—Rc,o)-P(Rc, O |R-c,o, ^o; t). 


(Bl) 


From the Smoluchowski equation, with the use 
of Eq. (1201) . it follows that in the mobility center, 
where D rt is symmetric, one has, 


dRc (Rc- 


Rc,o)^(Rc,B |Rc,o> ^o; t) 


n /* 

— J dRc (Rc - Rc,o)-P(R-c 1 ^ |R-c,o, ^o; t) 

= C rr j dR c (Rc - Rc,o)P(Rc, O |Rc,o, ^o; t) 
-2-^-D rt P(n\n 0 -,t), (B2) 

with P(fl\n 0 ;t) defined in Eq. (15U1) . In Eq. (IB2I) . 
the terms which contain ^ • D rt and D tr ■ ^ give 
identical contributions. 

Solving Eq. (IB2I) . we obtain 


= -2 / dr exp[C rr (t — T)}— — • D rt P 
. da 


(n 


^o;tJ 
(B3) 


We now insert the expression (IB3I) into Eq. dm, 
take into account that the Smoluchowski operator 
is self-adjoint, and perform the integration by parts 
with respect to a. We benefit from choosing the 
frame of reference, in which D rr is diagonal, with 


C rr n ( P ) = _yM u (p) ; (B4) 
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and we use the relations m, <EH) to write for a 
vector V and a matrix A, 

dV 

^4 = -(VxS) to . (B5) 

Finally, using Eq. (i38l) . we obtain Eq. (USD- 

Appendix C: Shift of the reference center 

The correlations {Au(t)Au(t )) 0 do not depend 
on the choice of the reference center. Indeed, they 
depend only on the rotational-rotational compo¬ 
nents of the diffusion tensor (or, equivalently, the 
mobility matrix). But these components do not 
depend on the choice of a reference center, as it 
follows from Eqs. (fTTl) - (fT3l) . 

In this appendix, we derive expressions which 
allow to transform the correlations of the rota¬ 
tional and translational Brownian displacements, 
(AM(i)AR(i)) 0 , from one reference center Ri to 
another, R 2 . We denote, 

R 2 = Ri -t- R 21 • (Cl) 

The difference R 2 i rotates with the particle, and 
therefore it can be expressed using the set of the 
orientation vectors u^ p \ with p = 1,2, 3, 

3 

R 21 = « (p) - (C2) 

P= 1 

Therefore, 

3 

AR 2 (t) = ARi (t) +J2 a P ( m(p) W - u(p) (0)) 

P= 1 

(C3) 

Since (Au(t))o = 0, the transformation of the 
rotational-translational correlations between the 
reference centers Ri and R 2 has the form, 

(Au(t)AK 2 (t)) 0 = (Aw(t)ARi(t)) 0 + 

3 

a v ( a u(t)u^ ( t . (C4) 

p =1 


To determine the difference, we take into account 
that Au(t) is proportional to see Eq. E51) . 

while according to Eq. (l26l) , u ^ (t) contains both 
TZ (a) and n {s \ 

u {p) (t) = 7 l {a \t) • u (p) (0) + n is \t) • tt (p) (0). (C5) 

However, as shown in Eq. the symmetric and 
antisymmetric parts of the rotation matrix are not 
correlated. Therefore, owing to Eq. (051) . the only 
contribution to u^ ( t ) in Eq. (IC4I) comes from 

K ia \t) • m ( p) ( 0) = A u(t) x «W(0). (C6) 

Therefore, 

(Au(t)u (p )(t)j o =(Au(t)Au(t)) 0 x u^( 0), (C7) 

and we finally obtain the transformation relation, 

(Au(t)AR 2 (t)) 0 = (Au(t)ARi(t)) 0 + 
(Au(t)Au(t)) 0 x R 2 i(0). (C8) 

This relation can be used to transform the ex¬ 
pressions dHU-dUD, valid for the mobility cen¬ 
ter Ri = Re, to account for the rotational- 
translational correlations determined for an arbi¬ 
trary center R 2 . 

The analogical transformation relations for 
the translational-translational correlations 
(AR(f)AR(f)) 0 contain also correlations in¬ 
volving ' KS s \t ), and therefore they cannot be 
expressed only in terms of the correlations C (t), 
and as such, they are not practical if only C (f) is 
measured. However, these additional correlations, 
(AR(f)u^(i)) 0 and (t)u^ (f)) Q , can be 

also determined in experiments. In this case, the 
corresponding expressions would be useful and 
can be easily derived based on the framework 
constructed in this work. 
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